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We studied the spin-dependent quantum transport properties using a simple modelling of a Datta- 
Das spin transistor. We refine previous results by accounting the propagation medium changes of 
opacity felt by itinerant electrons, when the gate-voltage is switched on and modelling them via 
the transversal energy levels mismatch. Monitoring the topological-dependent conductance, we 
are able to identify the device operating points. If the incoming electrons energy approaches the 
biased-induced barriers height, the spin-resolved conductance oscillations become significant. In a 
zero temperature picture, our computations of the spin-dependent conductance as function of the 
electric field at the region below the gate electrode suggest the feasibility of the modeled device. 
Although we demonstrate that phase time may not be spin-resolved, our simulation allows us to 
evaluate the time that takes an electron to experience a spin-flip process, resulting in an order of 
magnitude lower than typical values of the spin relaxation times. 


PACS numbers: 71.70.Ej, 85.75.Hh, 03.65.Xp, 85.75.-d 


I. INTRODUCTION 

One of the biggest challenges in Spintronics 0,0 is 
to achieve an efficient manipulation of the spin degree of 
freedom at zero mag netic field. Rashba spin-orbit cou¬ 
pling (SOC-R) 0 , 0 ] is among the most promising phe¬ 
nomena. SOC-R is a natural consequence of an asym¬ 
metry in the confining potential of carriers in low di¬ 
mensional semiconductor systems, so-called structure in¬ 
version asymmetry (SIA) [^. To intuitively understand 
SOC-R some authors consider that electrons feel a po¬ 
tential gradient due to SIA, which Lorentz transforms 
into an effective magnetic field that acts over theirs spin. 
Nonetheless, attention most be paid since Lorentz trans¬ 
formation neglects the atomic cores contribution to SOC- 
R, felts by a Bloch electron in a solid 0. 

The strength of SOC-R can be tuned by external elec¬ 
tric fields [ 3 , which imply a gate mechanism for spin- 
based transistors |7Hl^ Datta and Das spin field ef¬ 
fect transistor (SFET) Q (Fig. [1]) is a theoretical device, 
conceptually similar to an electro-optic modulator. Two 
ferromagnets act as polarizer and analyzer. The prop¬ 
agation medium between them —capable of inducing a 
gate-controllable net rotation of the spin orientation via 
SOC-R— is a two-dimensional electron gas (2DEG). 

Spin-dependent conductance in SFET-like systems has 
been widely addressed [TTI - [l8l| . In 2001, Mireles and Kir- 
czenow |ll| perform calculations of spin-dependent ballis¬ 
tic transport properties in quantum wires, in the presence 
of SOC-R. Using a tight-binding scheme, they showed 
that a strong SOC-R may lead to dramatic changes in 
the transmission of electrons. The transport of holes in 
p-doped hybrid (magnetic and nonmagnetic) structures 
was considered by Pala et al. [EH, exhibiting clear 
possibilities of current manipulation. The effects of the 
temperature [l^, external magnetic fields [T^ and the 
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direction of magnetization at the polarizer and analyzer 
electrodes jEllalEl have also been addressed. Recently, 
Gao et al. |16| carried out a simulation of a SFET, based 
on the nonequilibrium Green’s function formalism self- 
consistently coupled with a Poisson solver to produce the 
device I — V characteristics, obtaining good agreement 
with the recent experiments. However, these studies does 
not allow for the modifications in the opacity of the prop¬ 
agation medium that itinerant electrons feel when the 
gate-voltage is switched, as we will discuss later. 

Another issue insufficiently investigated is the spin- 
dependent tunneling time in SFET-like systems. As far 
as we know, only few studies have focused the quantum 
transport time under SOC-R in ferromagnetic/semicon¬ 
ductor/ferromagnetic systems. Based on the group veloc¬ 
ity concept, Wu et al. [1^ showed that as the strength of 
the SOC-R increases, the traversal time considerably de¬ 
creases and exhibits step-like behavior when the length of 
the channel growths. However, it should be pointed out 
that SOC-R appears in the system they are modelling 
as a consequence of an asymmetry introduced at the fer¬ 
romagnetic/semiconductor interface to fix some known 
problems related with the efficiency of the spin-injection. 
Then, as the potential gradient in this case is parallel to 
the transmission direction, the effective magnetic field is 
null. Hence, there is not spin-precession, in contrast to 
SFET-like systems where the potential gradient induced 
by the gate electrode is perpendicular to the transmission 
direction. 

This contribution focuses on the quantum transport 
properties in SFET-like devices. Using a simple mod¬ 
elling, combined with a suitable algorithm of solution, we 
show that spin-depend conductance exhibits the general 
features reported before, but strongly modulated by dis¬ 
persive effects, appearing when the SOC-R coupling pa¬ 
rameter is tuned via external electric fields. Although we 
demonstrate that the transport time —using the phase 
time concept— may not be spin-resolved, our simulation 
allows us to evaluate the spin-flip time, which is an im- 
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FIG. 1. (Color online) (a) Sketch of the Datta and Das SFET. 
Two ferromagnets act as analyzer and polarizer, the prop¬ 
agation medium —capable of inducing a gate-controllable 
net rotation of the spin orientation via SOC-R— is a two- 
dimensional electron gas (2DEG). (b) The five regions we 
are considering: Fl and Fa correspond to the ferromagnets, 
whereas 7, III and II stand for the regions free and the one 
right below the gate electrode, respectively. 


portant time scale for spintronics applications, together 
with the spin relaxation and decoherence times. 

The remaining part of this paper is organized as fol¬ 
lows: Sec. El is devoted to describe the physical sys¬ 
tem and the theoretical formalism describing electrons 
in SFET-like systems. The general properties of the sim¬ 
ulation procedure are presented in Sec. Uni In Sec. Iv] we 
discuss numerical results of spin-dependent conductance 
and phase time. Finally, in Sec. El we underline some 
concluding remarks. 


II. THEORETICAL MODEL 


We simulate the device depicted in Fig. [T] by a five 
regions system: Fl and Fr correspond to the ferromag¬ 
nets, whereas I {HI) and II stand for the free and gate- 
managed regions, respectively. The 2DEG is confined 
at the interface of a heterostructure, InAlAs/InGaAs for 
instance, where SIA lead to Rashba coupling. Spin- 
polarized electrons are injected by an ideal ferromagnet 
from the region Fl- Another ideal ferromagnet allows 
identifying the spin direction of the right side outgoing 
electrons. By “ideal” we overlook the actual limitations 
of ferromagnets [l^, [2^ , fixing the initial electrons spin 
polarization instead. For simplicity, we choose the same 
effective mass along the whole system. 

The [a;?/]-plane total Hamiltonian for a 2DEG confined 


on it, has the form 

+K) + Ebs ianK + oiRky 
-iank^ + aRky -;^{kl F k^) + Ebs 

( 1 ) 

Here we have introduced an extra term Ebs, standing 
for the bounded states at the 2 direction. In the infinite 
triangular well approximation to the heterostructure con¬ 
fining potential |2l|, [22| , Ebs are given by 


Eb 


s n 


= Cn 



1/3 


with: n = 1, 2, 3,..., (2) 


where Cn are zeros of Ai Airy function {Ci = 2.338, 

C 2 = 4.088, Ca = 5.521, ...), and A = ^ + Fext 
the electric field due to SIA. Note that E involves both 
intrinsic —density (n^) dependent— and external (Fext) 
contributions. 

The off-diagonal terms in ([T]) correspond to the SOC- 
R, where aR is the coupling parameter [Ej | 

ar,(F) — F ^ _ ( 3 ) 

> 2m* Eg {Eg + Aso){SEg + 2Aso)' ^ ’ 

Eg and Aso are the gap and splitt-off energies, respec¬ 
tively. 

Because of Rashba coupling, the eigenvalue prob¬ 
lem for the Hamiltonian Q H'P = E\F produces two 
branches in the dispersion relations 

= -^—k^+aRk, 

Zm* Zm* 

( 4 ) 


where k = Jk^ + ky, and eigenvectors 
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with 6 = arctanA:a;/fcj,. 

Hamiltonian ([T]) is time-reversal invariant, therefore it 
can not support spontaneous spin polarization of the elec¬ 
tron states [HI- However, it is capable of removing the 
spin degeneracy for k ^ 0 through an effective magnetic 
filed Bgff = fcy, fca;, 0), with fiR as the Bohr 

magneton. 


III. MULTICOMPONENT SCATTERING 
APPROACH 

Weakly coupled multi-mode quantum transport of 
electrons has been typically treated through one- 
dimensional single-mode approximations. The corre¬ 
sponding solutions, along with the envelope function ap¬ 
proximation, are multicomponent wave functions. How¬ 
ever, it is usual [l^, to arbitrarily cancel all com¬ 
ponents but one in order to compute transmission co¬ 
efficients, which may lead to lose valuable physical in¬ 
formation. Recently, an alternative formalism, named 
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Multicomponent Scatte ring Approach (MSA), has been 
proposed (See Refs. [2bl - l^ for details). MSA takes into 
account —jointly and simultaneously— all the propagat¬ 
ing modes, then a multicomponent description of trans¬ 
mission amplitudes is carried out in a natural way. 

Let 


2N 2N 

F{y) = CjFj(y) ( 6 ) 

j=i i=i 


be the general solution of the system described by the 
Hamiltonian ©• To introduce the transfer matrices we 
define the bi-vector 


' 5 '( 2 /) 


F{y) \ 
F'{y) ) 


( 7 ) 


and the state vector 


Hy) = 




( 8 ) 


fl' include the wave vectors and their derivative. The 
coefficients a and b are matrices (2 x 2 ). ^( 2 ) and V('^) 
are bidimensional vectors, whose components describe a 
propagating or evanescent mode, according to the energy. 

Based on the above definitions and taking into account 
that solutions (O are plane-waves, it is possible to estab¬ 
lish the crucial relation 


^{z)=M^{z). 


( 9 ) 


The matrix Af depends on the specific Hamiltonian, and 
is build up using the property (e““)^ = ae“^. In our 
particular case we have 


/I 0 1 0 \ 

I 0 1 0 1 I 

I ik 0 —ik 0 I ’ 

\o ik 0 —ikj 


( 10 ) 


being k = (l/h)\/2m*E. 

The first-kind transfer matrix Mfd {fd stands for func¬ 
tion and derivative), which relates the solutions and their 
derivative 0 between two points of the system and 
?/r, is defined by 

^{yn) = MfdiyR,y-L)^{yL)- (H) 

Likewise, the second-kind transfer matrix Msv {sv stands 
for state vector), which relates the sate vectors 0 , is 
defined by 


is obtained, which relates the two types of transfer ma¬ 
trices already defined. 

Following the transfer matrix Msv and scattering ma¬ 
trix S definitions we can establish 



respectively, where t and r (t' and r') are transmission 
and reflection amplitudes for left (right) incident par¬ 
ticles. Using (HI and ()15|1 together with the imposed 
boundary conditions it is possible to relate the scatter¬ 
ing amplitudes for the symplectic casefTzi [30l| 


t = a — f3S ^7 t' = d ^ 
r = r' = (35~^ 


(16) 


The connection between Transfer Matrix and Scatter¬ 
ing Matrix formalisms is the most remarkable feature of 
MSA. 

From relations (flHl) we could obtain relevant transport 
magnitudes. Considering the incidence of particles from 
the left only, the transmission coefficient from channel j 
to channel i is given by 

(17) 

Furthermore, they allow to compute the conductance 
through the channel i 




i=i 

the two-probe Landauer conductance 

y N N 

h 




i=l j=l 

the phase-transmission amplitudes 

Im tjj 


(j)ij = arctan 


Re ti 


(18) 


and the phase-transmission times or group delays 


(19) 


( 20 ) 


( 21 ) 


$(?/r) = Ms^(?/b,?/l)^(2/l). (12) 

These matrices satisfy the relevant physical properties 
and symmetries of the Hamiltonian (fT|). 

Using equations 0- 100 the transformation 

Msv{yR,yL) = Af~^Mfd{yK,yeW, (13) 


IV. SPIN-DEPENDENT TRANSPORT 
PROPERTIES 

Our gedanken experiment consists of introducing spin- 
polarized electrons in the three-region channel with SOC- 
R (I, II and III in Fig. [T] (b)), under normal incidence 
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FIG. 2. (Color online) SOC-R coupling parameter a_R @ 
(panel (a)) and gate-voltage induced potential barrier V (1221) 
(panel (b)) as function of the electric field at the region II, 
in units of the one at the bare regions (I and III). 


{kx = 0). We are able to tune up the coupling param¬ 
eter Or in the central region (region II, below the gate 
electrode), the length L of this region and the energy E 
of the incident electrons. We set the length of the re¬ 
gions / and III as Z = 10 nm. SOC-R in these regions is 
solely given by the intrinsic asymmetry of the confining 
potential, which depends on the surface charge density, 
TLs = 10^^ cm“^, from now on. Following the Landauer 
picture [3l|, we assume the whole system connected to 
infinite reservoirs of charge at different chemical poten¬ 
tial, avoiding the inclusion of an external electric field to 
move electrons along the system. 

There are two direct paths connecting electrons with 
the same spin {e^ and ej, e^) and two crossed 

paths connecting electrons with different spin (e-|- => ej, 
and ej, => e-j-). The availability of the last two is given 
by the SOC-R. In accordance with our consideration for 
regions F), and Fr, we are only interested in the paths 
ef ^ e-j- and Cf => ej, because we are injecting in the 
channel only spin-up polarized electrons. 

Taking as energy reference the first transversal state 
Ebs —given by Exp. m- electrons traveling from region 
I into II experience a potential barrier of height 



and thickness L, due to the mismatch of the energy lev¬ 
els between those regions, induced by the external elec¬ 
tric field. Then, the potential V essentially depends on 
the voltage at the gate electrode (region II, see Fig. 
[T]). Hence, by tunning Or the height of this barrier also 
changes, as shown Fig. [21 The coupling parameter 
growths following (see panel (a)), while V growths 
following (ESI) (see panel ( b)). 

It is a well-known result (lll - ll^ that spin-resolved con¬ 


FIG. 3. (Golor online) Closed (a) and Opened (b) SFET op¬ 
erating points: Fixing the polarization of incoming electrons 
as ^ = I 4)21 a spin-up G-f or spin-down polarized con¬ 
ductance is obtained at L = 600 nm and E = 0.48 eV for 
Fii = Fi III = ^ (panel (a) and Fn = 2.2Fi m = 2.2^ 
(panel (b)). 


ductance oscillates with the lengthening of the channel, 
due to the SOC-R-induced spin precession phenomenon. 
In fact, assuming that electrons move freely (fill no po¬ 
tential barrier (I22II ) along the positive y direction in the 
transport channel, following (0])-®, the equivalent wave 
function would be 

«’<!') = (!) + 71 (i.) (23) 


where kyi and ky 2 are the positive solutions of equations 
O- Consequently, the probability of detecting spin-up 
or spin-down electrons at a given L reads m 


F,p = |((10)|'®'(L))|2=cos2 


-{ky 2 kyl) 


(24) 


Fdo«,n = |((0 1)|^(F))|" =sin2 


2 iky2 kyl) 


(25) 


respectively. The oscillations are then straightforward 
from the presence of cosine and sine functions. This be¬ 
havior is observed in our calculations, as a general trend 
(see Fig. [3]). 

Both SFET operating points are shown in the pan¬ 
els (a) and (b) of Fig. |3l respectively. Fixing the po¬ 
larization of incoming electrons along the ^ direction 

^ = I for convenience, at null gate-voltage 

the intrinsic asymmetry of the heterostructure potential 
ensures that spin-up conductance (associated with 
the path e-j- e^) vanishes for L = 600 nm (Fig. |3| (a), 
red solid line). Hence, by conservation low, the spin- 
down polarized total outgoing flux maximizes ~ 1, 
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Fig. [3] (a), blue dashed line) but being associated with 
the path Cf ej,, becomes forbidden by the ferromag¬ 
netic analyzer, and then the device is closed. On the 
other hand, for a suitable value of the gate-voltage that 
strengthens the electric field 2.2 times respect to the one 
at the bare regions (F/,/// = the conductance for 

spin-rtp polarized electrons is magnihed (G-i- 1 and 

Gj, ~ 0), as shown in Fig. [3] (b) (red solid and blue 
dashed lines, respectively). In this case, the ferromag¬ 
netic analyzer allows the transport, thereby the device is 
opened. This two-states behavior of the SFET had been 
proposed for codifying and manipulating information 
(“ 0 ” and “ 1 ” values of a bit) instead of the current field 
effect transistors. 

At a given value of the Fermi energy, by evaluating kyi 
and ky 2 in Exps. ([4]) we have Ei{kyi) = i? 2 (fcy 2 ), then 


2 m* 




E2{ky2) — Ei{kyi) — 0, 


fcyl) + aR{ky2 + kyl) = 0, 


then ky 2 


kyl — 


2m*aR 


(26) 


From (I^Sl) it is straightforward to note that Pup (EH) and 
Pdown (l25l) do not depend on the energy [Tl|, therefore 
the consequences of changing the energy on the spin- 
resolved conductance are typically leaved out. In our 
modelling, the energy plays a crucial role because we do 
consider the dispersive effects of the constriction at the 
region II (see (l 22 |l and its discussion). Actually, the 
additional “noise” —fast oscillations— in spin-resolved 
conductance curves, at finite values of the gate-voltage 
(Fig. [3] (b)), is clearly related with this. To evaluate 
the significance of those effects, in Fig. [4] we plot the 
spin-resolved conductance in a particular range of en¬ 
ergy, setting all parameters as in Fig. [3] (b). The gate- 
voltage induced potential barrier (( 22 l) has the height 
V{Pii)+Ei,s = 0.37eV, which is lower than the whole en¬ 
ergy range considered, then in principle is not tunneling 
but transmission the quantum phenomenon we are con¬ 
sidered. Note that the total conductance, G.f+ ^ 4.1 could 
be reduced from 75% down to 25%. The existence of res¬ 
onant peaks is a well-known phenomenon, associated to 
the condition k{E)L = mr, where n is an integer. 

Fig. [5] displays the spin-resolved conductance as func¬ 
tion of the electric field at the region II (induced by the 
gate electrode), in units of the one at the bare regions (/ 
and III). This [I— E)-like curve characterizes the behav¬ 
ior of the modeled SFET. Note that for F 77 = Fjjji = 
e^nsjene the outgoing flux is spin-down polarized, leading 
to a closed state, while for F 77 = 2 .2Fijii = 2 . 2 e^ns/eoe 
is spin-itp polarized, leading to an opened state. We chose 
the energy of the incoming spin-np polarized electrons 
as E = 0.48 eV, in order to avoid the zone of large- 
amplitude oscillations of the spin-resolved conductance 
(see Eig. |3|). Even that, for the highest values of the 
electric field (F77 > 1.6^7 777) those oscillations are ob¬ 
served. 



0.38 0.40 0.42 0.44 0.46 0.48 0.50 

E [eV] 

FIG. 4. (Color online) Spin-resolved conductance as function 
of the energy for the conditions represented in Fig. [3] (b), at 
L — 600 nm, Fu = 2.2Fijii and V{Fji) E^s = 0.37 eV. 



FIG. 5. (Color online) Spin-resolved conductance as function 
of the electric field at the region II (induced by the gate elec¬ 
trode), in units of the one at the bare regions (/ and 111). 
L = 600 nm and E = 0.48 eV 


We are interested now in evaluating the quantum 
transport time in the system we are modelling (see 
Eig. HI). To answer how long it takes an electron to pass 
through a quantum system is still a controversial question 
[ 3 ^ . Phase time, defined as in (I3T]1 . is a well-established 
definition of this time, and has been extensively used in 
researches regarding quantum tunneling [271 - 1^ [3413^ . 

The first task then is to compute the phase as func¬ 
tion of the energy following (1301) . Although not shown 
here, if one plots the phase for the allowed paths, (j)^ 
(e^ ^ e^i) and (^ 4 , (e^- => ej,), one may note that the 
corresponding curves are shifted in the energy. This is 
consistent with the fact that electrons with different spin 
polarization propagate along the system with different 
quasi-momentum, as in the “toy model” described by 
the solution (l23l) . Consequently, those curves have equal 
derivative and the phase time is the same along the direct 
and crossed paths. This may be explained taking into 
account that manipulation of spins in a SEET via SOC- 
R —in itself— is not a dispersion-selective phenomenon 
thus the coherence is preserved. The opposite picture 
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FIG. 6 . (Color online) Phase time as function of the electric 
field at the region 11 for two particular values of the energy, 
E — 0.40 eV (squares) and E — 0.48 eV (circles). L = 
600 nm. 



FIG. 7. (Color online) Phase time as function of the length L 
of the region under SOC-R (region II) for both closed {Eu = 
Eiji, circles) and opened {{Fu — 2.2Fiji, squares)) states, 
with E — 0.48 eV. 


occurs in ferromagnets or dilute magnetic semiconduc¬ 
tors, where spin separation features in time have been 
reported [H, A simple modelling of those systems 
[la [T^ predicts that they behave as a potential barrier 
for a particular spin polarization and as a potential well 
for the other, resembling an effective Zeeman splitting. 

Fig. [6] displays the phase time r as function of the elec¬ 
tric field at the region II for two particular values of the 
energy, E = 0.40 eV (squares) and E = 0.48 eV (cir¬ 
cles). At first sight, one may note that the phase time 
is lower for the higher energy and its value increases as 
a general trend in both cases, when one lets the electric 
field growths. This is intuitively understandable noting 
that the height of the potential barrier (1221) also growths, 
as shown in Fig. [H (b). However, no significant changes 
are observed in r for £1 = 0.48 eV when one switches the 
device from the closed state (F// = Ejjn, depicted in 
Fig.|3](a)) to the opened state {Eu = 2.2F/j//, depicted 
in Fig. [3] (b)). On the other hand, for E = 0.40 eV, r 
varies from ~ 20 ps to ~ 45 ps, which may be quite un¬ 
desirable for practical applications. Another interesting 
issue here is the presence of oscillations in both curves, 
being more evident for E = 0.40 eV. These oscillations 
have been reported before in the case of holes [H, |2^ , 
but as a function of the energy for values above the po¬ 
tential barrier, therefore they have the same nature that 
the ones we are reporting here. 

In Fig. [7] we plot r as function of the length L of 
the region under SOC-R (region II) for both closed 
{Fjj = Ejjj, circles) and opened {Eu = 2.2F/j/, 
squares) states, with E = 0.48 eV. Note that for L = 
600 nm, which is a kind of spin-flip length for Ug = 
10“^^ cm~^ (see Fig. [3] (a)), the phase time for the 
opened and closed states barely differs in ^ 4 ps. While 
the phase time growths linearly with L for Eu = Ej u 
(zero gate-voltage), exhibiting a classical-like behavior 
{tf = L^Jm*/{2E)), it follow a non-trivial increment 
when Eu = 2.2Ejjj. From Figs. [S] and [7] we may infer 


device we are modelling are ~ [10 — 50] ps, while the spin 
relaxation time is typically in the order of nanoseconds 

H. 

It should be mentioned that is the first time that MSA 
formalism deal with systems of hundreds of nanometers, 
which is the typical scale length where spin-flip occurs 
[HI [alii- So far, the magnitudes (fT71) - (l^ had been 
successfully studied using MSA in systems with length 
ranging form ^ 10 A up to ~ 100 A. General properties, 
as the probability flux conservation through our system, 
was verified during the simulations. All physical param¬ 
eters not mentioned above are referred to Ino. 53 Gao. 47 As 


V. CONCLUSIONS 


The changes in the Datta-Das spin-transistor propaga¬ 
tion medium opacity, given by the manipulation of the 
SOC-R strength, may considerably reduce the conduc¬ 
tance of the device. The spin-resolved conductance have 
exhibited oscillations as expected, but modulated by dis¬ 
persive effects when the gate-voltage was turned on, lead¬ 
ing to significant reduction of the total tunneling con¬ 
ductance. Nonetheless, in a zero temperature picture, 
present simulation allows us to identify suitable topo¬ 
logical and external parameters, to obtain the operating 
points of the SFET and to prove its feasibility as sub¬ 
stitute for current field effect transistors. Although we 
demonstrate that phase time may not be spin-resolved, 
we were able to evaluate the time that takes a spin-up po¬ 
larized electron to switch it down, resulting in a range of 
~ [10 — 50] ps, an order of magnitude lower than typical 
values of the spin relaxation times. The systematical pro¬ 
cedure carried out here, provides the basis for studying 
novel SFET-like configurations 0 and/or addressing p- 
doped SFET-like systems where the phenomenology is 
strikingly different respect to the electronic case [^. A 
work devoted to those issues is in progress and will be 
published elsewhere. 
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